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Abstract 

The transition pion-photon form factor is studied within the framework of Light-Cone QCD 
Sum Rules. The spectral density for the next-to-leading order corrections is calculated for any 
Gegenbauer harmonic. At the level of the next-to-next-to-leading (NNLO) radiative corrections, 
only that part of the hard-scattering amplitude is included that is proportional to the /3-function, 
taking into account the leading zeroth-order harmonic. The relative size of the NNLO contribution 
in the prediction for the form factor ^^{Q'^) has been analyzed, making use of the BLM scale- 
setting procedure. In addition, predictions for the form factor F'^ are obtained that turn out 
to be sensitive to the endpoint behavior of the pion distribution amplitude, thus providing in 
connection with experimental data an additional adjudicator for the pion distribution amplitude. 
In a note added, we comment on the preliminary high-Q^ BaBar data on F'^ '^'^ arguing that the 
significant growth of the form factor between 10 and 40 GeV^ cannot be explained in terms of 
higher-order perturbative corrections at the NNLO. 
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I. Introduction 



Although higher-order calculations in QCD perturbation theory have already a long his- 
tory, httle is known about exclusive processes at the next-to-leading order (NLO) level 
B i, i i, i i, and beyond 0, i, because these are quite complex in detail. In view of 
more and more high-precision experimental data for a variety of hadronic processes becom- 
ing gradually available, the importance of such higher-order calculations exceeds the pure 
theoretical interest and acquires phenomenological relevance. In particular, processes with 
two photons in the initial state, one far off-shell and the other quasi real, 



7 +7 



provide a useful tool to access (after their fusion) the partonic structure of the produced 
hadronic states, e.g., pseudoscalar mesons. 

Experimentally, the photon-to-pion transition form factor within this class of two-photon 
processes has been measured by the CLEO Collaboration [9| with high precision and ex- 
tending the range of up to 9 GeV^, as compared to the previous low- momentum CELLO 



data [10|. Theoretically, this high precision allows one to test models and fundamental 
quantities, like the pion distribution amplitude (DA), the applicability of QCD factoriza- 

and 



tion, etc. — see jl, 0, 1111 
references cited therein. Moreover, one can determine [26| a compatibility region between 
the CLEO data and constraints derived from lattice simulations on the second moment of 
the pion DA 28|, |29|. This information can then be used to extract a range of values of the 




fourth moment of the pion DA that would simultaneously fulfil both constraints (CLEO and 
lattice). This prediction j26| can provide a guide for the determination of this moment on 
the lattice, a task that has not been accomplished yet. 

For two highly virtual photons, perturbative QCD works well because factorization at 
some factorization scale /i| applies, so that the process can be cast into the form of a 
convolution 



2 2 \ 



(x,/x^)+0(Q- 



which contains a hard part C, calculable within perturbation theory, and a wave-function 
part (y9,r that is the (leading) twist-two pion distribution amplitude |30| and has to be mod- 
eled within some nonperturbative framework (or be extracted from experiment). Here, the 
omitted twist-four contribution represents subleading terms in the operator product expan- 
sion (OPE), which are suppressed by inverse powers of the photon virtualities. 

To be more precise, consider the hard process of two colliding photons producing a single 
pion, 7*(q'i)7*(q'2) ^°(p)) which is defined by the following matrix element [17 



d xe~ 



■iqi-z 



(TT^ip) I T{j^{z)UO)} I 0) = ^e^,„^grg2^F^'^^-(Q^g 



:i.2) 



where 
the momentum p 



= — ^2 denote the virtualities of the photons, vr'^(p) is the pion state with 
Qi + q2, and = {^ifj^u — ^d'j^d) is the quark electromagnetic current. 
This process is illustrated graphically in the left panel of Fig. [T] and has been examined 
theoretically, for instance, in 

If both virtualities, and g^, are sufficiently large, the T-product of the currents can 



be expanded near the light cone 



0) by virtue of the OPE to obtain the well-known 



2 



pQCD 



LCSR 



FIG. 1: Pion-photon transition form factor in pcrturbative QCD (left) and using Light Cone Sum Rules 
(LCSR) (right). The right graphics illustrates the situation when one of the two photons is real and 
perturbation theory becomes inapplicable, because the hadronic content of the photon starts to be relevant. 



leading-order expression for the convolution in Eq. (II. ip 



Jq Q^x + q^x 



:i.3) 



where we have used the abbreviations x = \ — x and 



V2 



:i.4) 



In contrast, the kinematics probed in the CELLO [lo| and the CLEO [9| experiments in- 
volves a quasi-real photon with — > 0. At such a low virtuality, the hadronic content of the 
quasi- real photon, i.e., its long-distance structure 3] becomes important (see the diagram 



in the right panel of Fig. [T]), thus preventing a straightforward QCD calculation [13| of the 



form factor F"'*"''^ [Q"^ ^ 0) = F'^*'^'^ [Q'^) on the ground of factorization. The method of 
Light-Cone QCD Sum rules (LCSR for short) allows one to avoid this problem by providing 
the means of performing all QCD calculations at sufficiently large (7*) and then use a 
dispersion relation to "approach" the mass-shell photon (7) with zero virtuality. This cal- 
culational scheme, which can accommodate the large-distance properties of the photon, i.e., 
its hadronic content, was proposed by Khodjamirian in [3] and the form factor F^*^'^{Q'^) 
was calculated at the LO level of the LCSR including also twist-four contributions. 

The core ingredient of the LCSR approach is the spectral density, which provides a 
powerful tool for a quantitative description of hadronic processes in QCD in terms of a 
dispersion relation: 

00 
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77r 



ds 



pp'^^"(g' ,g) 

s + q^ 



:i.5) 



[Note that here the label "phen" abbreviates phenomenological] . Effects due to the long- 
distance dynamics of the 7*7 — > vr" process are partly contained in the form factor F'^ '''^{Q'^) 
(which can be obtained by means of quark-hadron duality) and also in the tt distribution 



amplitudes of different twists [3l|, |32|, |33|, |3J] . Within the LCSR approach the form factor 
p"i ptt^q2^ appears inevitably because one assumes that the spectral density, entering the 



dispersion relation, can be approximated by the ground states of vector mesons [ij, l35|, like 



3 



the p and the uj} 

At this point, two important remarks are in order, (i) Contrary to previous calculations, 
e.g., in jl, 6|, we will not adopt a zero- width approximation here, but use instead a more 
realistic Breit-Wigner ansatz for the effective resonances (see Sec. |V]). (ii) The scaled form 
factor Q^F"* '"^(Q^), obtained in the framework of LCSRs, depends at large mainly on the 

^^^{x)\^=,, with e ~ ^ 

around the origin, where sq is the duality interval entering the model for pP^™. This feature 
appears to be opposite to the case of the Q'^F^*^^{Q'^) (scaled) form factor, that depends 
mainly on (though it is not directly proportional to) the inverse moment 221] (x"^),r = 
V^7r(a;; fi^)x~^dx, cf. Eq. (11. 3p evaluated at the scale — > 0, because the latter is an integral 



differential pion characteristic —LpT,{x)\x=e with e ~ — <^ 1 being a small neighborhood 



characteristic of the pion DA [18|, [l9[. Hence, Q^F^*^p'^{Q'^) can provide complementary 
information on the pion DA and help discriminate among various proposed pion DA models. 

The structure of the paper is the following. In the next section, we recall the formalism 
of LCSRs for the form factors F"'''"''^{Q'^), F'^*p^{Q^) and construct the spectral density in 
a systematic way. This calculation is extended beyond the LO in Sec. IIIII and an explicit 
expression for the spectral density at the NLO for any index n — the latter indicating the 
order of the expansion in Gegenbauer harmonics — is derived. The further extension to the 
NNLO is also given in this section. Actually, we include only the /^o-proportional contribu- 
tions that can be obtained from the corresponding terms of the hard-scattering amplitudes 
and denote it by NNLO^. The effects of the NNLO contributions are discussed in Sec. IIVI 
in connection with the BLM prescription (and its modifications js^). Our predictions for 
the form factors F'^ ^^{Q"^), -F'^*'"^(Q^) are presented in Sec. |Vl where we also provide a 
comparison of F^''^'^{Q'^) with the experimental data. Section IVl] contains our conclusions 
emphasizing our main results. Important technical details are provided in three dedicated 
appendices. In a Note Added, we point out that the new BaBar data on F^'^'^ , which show 
a significant growth of the form factor beyond 10 GeV^, cannot be described within the 
QCD convolution approach. 



II. FT^*7^(q2) i77VT(g2^ j]vj LIGHT-CONE SUM RULES. FORMALISM 

Here we present the theoretical description of the transition form factor F^'^'^iQ/^) of the 
exclusive process 7*7 — vr , employing the framework of light-cone sum rules [IJ, 53, Hi^ 



beyond the next-to-leading-order (NLO) of perturbative QCD.^ An integral part of this sort 
of calculation is the form factor F^*^'^(Q^), describing the transition 7*p — ^ vr, which will, 
therefore, be computed in parallel. 



^ For the sake of simplicity, one sets the masses of the two vector mesons equal and appeals to isospin 

symmetry to treat both particles in terms of a combined effective resonance. 
^ For the sake of clarity, we use the following notation for the form factors: The associated reaction is 

denoted by superscripts, whereas the calculational context, e.g., QCD, is marked by subscripts. Note that 

our notation differs from abbreviated notations used in the cited works. 



4 



A. Factorization 



The calculation of F'^'*"''^ [Q'^) proceeds through the following main steps: (i) First, the 
form factor FQ^^Yi^{Q'^, q^) at large Euclidean virtualities of the photons, Q^, > 1 GeV^ is 
calculated, (ii) Then, an appropriate realistic model for the spectral density at low s, based, 
for instance, on quark-hadron duality, is constructed, (iii) Finally, the dispersion relation 
for the form factor F"'"'* [Q"^ , q^) is exploited. 



Applying the factorization theorems [3l|, [32, l33|, l3J|, the form factor Fq(^p''((5^, q^) can 
be cast in the form 

^QCD''(<5^ (f') = ®ip^ + higher-twist contributions . (2.1) 
The hard-scattering amplitude T for this process, written below in the square brackets, 

F^ci'^iQ', = [To(g^ q'; x)+ a] T,{Q\ q'; fj^l; x) 

+ T2{Q^,q^;iJ,l;iJ,l;x) 

+ ■■■ ]^^^^\x;fil) 

+ higher-twist contributions, (2-2) 

a. = ^4^, ®^Cdx, (2.3) 

is calculable within perturbative QCD. The symbols /ir and /ip denote, respectively, the scale 
of the renormalization of the theory and the factorization scale of the process. The pion 
DA Lp^T^\x] fip) of twist two, entering the convolution with the hard-scattering amplitude, is 
inaccessible to perturbative QCD and demands the application of nonperturbative methods 
(see Section rVAl) . We quote the well-known result for FqJ^^ at LO in that also includes 
the twist-four contribution llJl, viz.. 



F^o'QCBiQ^l') = iVT{To(Q^g^;x)®^(2)(a;)- [To(Q^g2;x)]'®^(^)(x)| (2.4) 

with 

ToiQ', q') = 2 + ^ ^) (2.5) 

2 Q'^x + q'^x 

and where ip^^\x) is a naturally appearing combination of twist-four pion DAs (for more 
details, see j4o| and Sec. IV Ap . An explicit expression for Ti has been obtained in 0, 0]. 
More recently 01, the general structure of T2 was investigated and its /3-part, bo ■ Tf3, was 
calculated. For convenience, both amplitudes are included in Appendix [XI In fact, the 
hard-scattering amplitudes To,Ti,T2 also determine the spectral densities p^°\ p^^\ p^'^\ as 
one can see from the following generic expression 

Tm 

p{Q\ s) = — [{T ® {Q\ -s)] , s>0. (2.6) 

71 

The core issues of the hard-scattering amplitudes are listed below, while characteristic Feyn- 
man graphs for each order of the perturbative expansion are depicted for illustration in Fig. 

m 
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(a) 




LO y*^ NLO 



(b) 




(c) 



FIG. 2: Typical Fcynman graphs contributing at different orders of the perturbative expansion of the pion- 
photon transition form factor in QCD. (a) LO, (b) NLO, and (c) NNLO, with the momenta of the various 
particles being indicated explicitly in (a). Both colliding photons arc considered to be highly virtual. 



1. LO perturbative QCD: Tq, see Fig. |2](a). 

In this order of the expansion we have To((5^,g^; 
scale cannot be determined uniquely. 



x) ® ip{x; fip) and the factorization 



2. NLO perturbative QCD: Ti. A typical diagram is shown in Fig. [2](b). 

Here the hard-scattering amplitude starts depending upon the factorization scale /xp- 
In explicit terms this reads 

Ti {x; Q\ g2) ®^{x)= To{Q\ g^; y) ® {CfT^) (y, x) + Ln(y)y ^ (y, x)] ® ^ {x) 
and the /ip-dependence enters via the logarithm 



Ln(y)^ln[(g2y + g2y) 



(2.7) 



whereas the LO Efremov-Radyushkin-Brodsky-Lepage (ERBL) kernel 31, 13, 4l| V^^^ 
is given in the next section with more details provided in Appendix [Bl As regards the 
term T^^\y^x)^ defined in flA.7l) . it contains contributions from the partial kernel V'^ 
(which enters V'^^^ and is generated by contracting the diagram in Fig. [2](b)) and also 
the kernel g stemming from V'^^'> and defined in flA.5l) . 

NNLO perturbative QCD: T2 with an example depicted graphically in Fig. [2](c). 
In this case, things are more complicated. First of all, we reiterate that only the 
/5o-proportional contributions to T2 (termed Tp) will be considered, making use of the 
calculations performed in 0]. This is, because that part can provide the right size 
of the whole NNLO contribution using the scheme-independent ELM scale-setting 
procedure to eliminate the appearance of the /3-function in the perturbative series. To 
discuss the structure of T^, let us first present it explicitly (postponing details to the 
appropriate sections to follow): 



Tp®^= In Ti ® + To g 

-^Ln2(|/)yw| 



r t(2) 



CFLn(y) 



^13 



(2. 



Second, as one sees from this expression, the hard-scattering amplitude depends ex- 
plicitly also on the renormalization scale, /x^, owing to the renormalization of the 
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strong running coupling. As a result, additional logarithms of the form \n{fip/fi^) 
appear. These terms are controlled by the renormalizat ion-group (RG) equation and 
can be resummed by applying the BLM procedure 42| which amounts to a rescaling of 
the argument of the running coupling to another value. This will be discussed later in 
more detail. Third, the logarithms Ln(?/) and Ln^(?/), which bear the /Xp- dependence 
[cf. Eq. fl2.7p ]. are accompanied by elements of the kernel V^^^ and the kernel V^^^ 
governing the NLO and LO evolution, respectively. 



B. Construction of the spectral density 

We continue here with the systematic construction of the spectral density. In doing so, 
we will make use of another set of variables, (x, Q^), instead of the usual set (s, Q"^). In the 
LO approximation, p^^\x) follows from the definition (12.61) upon inserting in Eq. (12. 5p the 
explicit expression for To(Q^,g^;a;) [l, Then, we obtain 



Im 



TT 



To(Q',-s) +tw-4 



Q 



(2.9a) 
(2.9b) 



[cf. first item in Appendix [A]. One notices that p'-"^ is directly proportional to the pion DA of 
leading twist two, (fin \ and the derivatives with respect to x of the twist-four contribution, 
Lp^^\ This observation allows one to simplify the subsequent analysis by introducing the 
normalized spectral density p{x,Q'^) via Eq. fl2.9bl) . 

Expanding if^\x; fi^) in terms of the eigenfunctions ipnix) of the one-loop ERBL evolu- 
tion equation, which coincide with the Gegenbauer polynomials Cn^^(.^), we find 

ip^^\x;p^) = M^)+ Yl «n(/i')V^n(x) , V^„(x) = 6x(l-x)C^/2(2x-l). (2.10) 

n=2,4,... 

In this representation, all scale dependence is contained in the coefficients anip"^) and is 
controlled by the ERBL equation. For this reason, it is convenient to project the leading- 
twist part of the spectral density p(Q^,s) on the same basis of eigenfunctions {ipn} and 
introduce the partial density pn, which has the form 



Pn{Q^,x) 



where x = Q"^/ (Q^ + s) and 



Pnix) = (x) + a] pW(x) + al pf )(x) + . . . . 
Therefore, from definition (12.111) and Eq. (I2.9ap it follows 



(2.11) 



(2.12) 



p(°)(x)=V^„(x). 



(2.13) 
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C. Dispersion relation for F'^' 



■yiT 



The spectral density, discussed above, allows us now to construct the phenomenological 
spectral density 



(2.14) 



which consists of two parts. The first term, p^, (with "h" denoting hadronic) encodes the 
hadronic content of the spectral density, reexpressed in terms of the 7*p n transition form 
factor FT*''^: 



m 



(2.15) 



where we assumed that the p and u resonances have the same mass and can be represented 
by a (5-function. Later on, we are going to show how this simple ansatz can be improved to 
include a finite width IJ]. The second term in Eq. (12.141) represents its perturbative part 



n 



-^QCD (Q ' ^) 



(2.16) 



in which Fq^-^^ {Q"^ , q"^) will be computed according to Eq. (12.21) (i.e., in convolution form) 
including contributions up to the NNLO/j. 

Then, substituting pP^™ in the dispersion relation for F'^*'^'^((5^, q^) in (II. 5p 14| . we obtain 



F^*^-{Q\ q') = ^ff ^J '"f + 



ds 



(2.17) 



where we have assumed that the hadronic spectral density, p^, in the dispersion integral can 
be approximated by the expression 



ml + g2 



so 



s + g2 



ds , 



(2.18) 



which determines the structure of F'^* p'" (Q"^) . Here, and below, sq = 1-5 GeV^ denotes the 
duality interval in the p- meson channel. 

In order to derive the LCSR for F'^*p'^{Q^, M^) [l^, we first perform a Borel transforma- 
tion of Eq. (I2.18P with respect to the virtuality q^ and then insert the result into Eq. (12.171) . 
Finally, taking the limit ^ 0, we arrive at 



so 



ds Im 



-^QCD 



+ 



ds Im 

S IT 



-^QCD (Q ' ^) 



(2.19) 



where ^0.7 GeV^ is the typical value of the Borel mass parameter and rrip is the p-meson 
mass. The first term in Eq. (I2.19p . which is proportional to F'^'^'", expresses the hadronic 
content of the quasi-real photon, whereas the second term encodes the pointlike subprocesses 
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governed by QCD perturbation theory. The spectral density, given by Eq. (12.61) . will allow 

77*71" / 
lcsrI 



us to obtain both parts of -^lcsr(Q^)' 



F^clW) = Nt \ ^V(Q^ M') + ^,H{Q') } , (2.20) 



where 



so 

V{Q\ M^) = jp{Q\ s)e("'""'V''"ds , (2.21) 



and where 



Hm=Q'' ip{Q\s)- (2.22) 

in correspondence with our remarks below Eq. (12.191) . Employing this notation, F'^ '"^(Q^) 
from Eq. (l2.18p can be recast in the form 

F^''''^{Q'') = ^V{Q\M^). (2.23) 

For the sake of completeness and for future use, we also display the spectral images of Eqs. 

(I2.20p . (12.211) . and (I2.22p in terms of the normalized spectral density, defined in Eq. (]2.9bp : 

^(Q^M^)^ L(g-|i5)p(..«^)f, (2.24) 



Xo 



H{Q') = Jp{x,Q')-, (2.25) 


where Xq = Q'^/{Q'^ + Sq) and s = Q'^x/x. The spectral density pn beyond the LO will 
be constructed and discussed in the next section. Let us remark at this point that in the 
following we are going to consider also mixed forms of these sets of variables, i.e., a spectral 
density which depends on both variables s and x. This should not cause any confusion 
because it is understood that one has to replace each time the appropriate variable, i.e., 
either s = Q^xjx ox x = Q"^ / [Q"^ + s). 



III. STRUCTURE OF THE SPECTRAL DENSITY BEYOND LO 

This section extends our calculation of the spectral density beyond the leading order of 
perturbation theory. In the following two subsections we will take into account the radiative 
corrections at the NLO and also at the NNLO; in the latter case only the /3o-proportional 
terms will be included. 
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A. Spectral density in NLO 



Let us first write down the final result for for any index n and tlien proceed with a 
systematic discussion of its structure: 



1 



(x; s/f,l) = (-3 [1 + v\n)] + y - In' (|) + 2vin) In 



-2v(n)- 



du 



du 



3 \x 



U — X 



+ (x 



U — X 



In ( 1 - - 

u 



X] 



\x 



XI 



tpnix) (3.1a) 
(3.1b) 
(3.1c) 



The partial cases n = 0,2,4 coincide, after some algebraic manipulations, with the results 
obtained before by Schmedding and Yakovlev (SY) The above expression — besides 

being valid for any n — also reveals how the radiative corrections manifest themselves in its 
various terms, as we will now explain. 

The quantities v^{n) and v{n) in fl3.1al) and (13.1bp are the eigenvalues of the corresponding 
parts of the one-loop ERBL evolution kernel V^^\ notably, V^, and (V^^ + V^). These are 
given by 



V^'\x,y) = CF2[V^{x,y)+Vlix,y)] = C^V+ix^y) 
V^''\x,y)^My) = CF2v{n)Mx), 



where the eigenvalues of the partial kernels V"'' are obtained from 



1 



{n + l){n + 2) 



1 

2 



Vl{x,y)0My) = v''in)ij^{x), v\n) = 2[^{2) - ^{2 + n)] 



(3.2) 
(3.3) 



(3.4a) 
(3.4b) 



Here tpiz) = d\nT{z)/dz and V+(x,?/) = V{x,y) — 6{x — u) V{t,u) dt. The latter defi- 
nition reflects the vector-current conservation, while further details pertaining to the above 
equations are relegated to Appendix [Bl Note that the term (13.1bp originates from the log- 
arithm in Ti [cf. Eq. f l2.7p ]. Therefore, its contribution to the discontinuity contains, in 
comparison with p^^^ in fl2.9ap . a new element which is discussed in Appendix [Xj item 2, 
Eq. flA.121) . The corresponding kernel in Eq. fl3.1b[) maps each monomial term again onto a 
monomial one, e.g.. 



du 



u 



X 



u — X 



9{x > u) 



u — X 



On' 



[ilj{n 



^(l)]x" 



(3.5) 



where a notation deviating from the usual + prescription 



[f{x,u)]^^^^ = f{x,u)-5{x-u) / f{x,t)dt 



(3.6) 



The corresponding explicit expressions, denoted Aq, A2, and A4, can be found there. 
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has been used. Keeping in mind that, ultimately, we have to integrate the spectral density 
p{x) over X, it is particularly useful to recast Eq. (]3.1bp in terms of an expansion over an 
orthogonal polynomial basis, e.g., over the eigenfunctions {ipn}'- 



du 



U — X 



\X 



X] 



1=0,2,. 



3 
2 



(3.7) 



The expansion coefficients are given in Appendix A, item 2. 

As regards the term fl3.1cp . it originates from the kernel g{y,u), introduced in 43|] and 
discussed in (with some technical details being provided for the reader's convenience in 
Appendix [Al item 2). This kernel, termed g in 43|, is not diagonal with respect to the 
{'?/'„ }-basis and is responsible for the apparent breaking of the conformal symmetry in the 
MS-scheme j43|. Notice that also the kernel in Eq. (13. Id) maps each term ipn onto a sum of 
i^iiy), I < Therefore, (13. Id) can be cast in the form of the following algebraic expansion 



1 r 



1pn{u) - 1pn{x) 



u 



X 



In 1 



du 



[X 



X 



^ Gnl ^l{x) 



(3.8) 



1=0,2,. 



Note that the various terms in the spectral density p*-^^ have a one-to-one correspondence to 
the kernel V^^'^ and its elements, and partly also to the kernel V^^"* via its element g. 

(i) For the special case n = 0, iI)q{x) = Qxx so that the dependence of Po^^ on the 
factorization scale p| disappears, owing to the fact that the asymptotic DA does not evolve 
in this approximation. Indeed, in this case, the following chain of evident simplifications is 
induced 



^;''(0) = i;"(0) = v{0) = 0; 



u 



X 



u 



tpo{u) + {x ^ x) = tpoix) . 



(3.9) 
(3.10) 



-ix) 



Substituting Eqs. (13. 9p . (13.101) in expressions (I3.1al) - (l3.1cp . one arrives at 

r2 



-3-2 



In^ 



i^o{x) 



(3.11) 



This expression agrees with the result obtained in [sl for p^q\q'^, s). 

(ii) For the general case of an arbitrary n, the "nondiagonal" (in ^n) part of p^^^ in the 
second line of the expression (note that s = Q'^x/x) 



In 



{x-s,pI) = |-3[l-^;'^(n)] 

Xlpnix) - 2 



TT 

y 



In^ 



2v(n) In 



2v(n) In 



n 

E 

1=0,2,.. 



Gnl i^i{x) + v{n) ^ bnii^i{x) 



=0,2,. 



(3.12) 



has been rewritten in terms of the known coefficients Gni and bni, supplied in Appendix 
IX] in terms of Eqs. flA.91) and flA.13p . This way, we have achieved that Eq. (I3.12p is a 
purely algebraic expression for p. To return to the n = case, one should set in Eq. (I3.12p 
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v^{0) = v{0) = and Goo = 1- Expressions (13.1ap - fl3.1cp . or, equivalently, (13.121) . provide 
an effective tool for analyzing any model of the pion DA within the LCSR approach. At this 
point it is worth comparing the NLO contribution, originating from two different approaches, 
with respect to the case when one of the photons becomes real. First, we have the expression 
obtained directly from the factorization formula in Eq. (1A.7P for = 0, 0, 3|, i.e.. 



X 



X 



ln(x) - 9 



i^oix) 



dx 



X 



-15Cf . (3.13) 



Second, starting from the dispersion relation given by Eq. ( ]2.25p for H, one can return to 
the previous expression by setting sq = 0, Xq = 1 to get 



X 



dx = Cf 



^ - in^ r- 

3 \x 



X 



-15Cf 



(3.14) 



where we have used p^^^ from Eq. (13. lip . The outcome of both expressions is the same, 
though the integrands are different. However, the In^-terms in both formulas (I3.13P and 



(I3.14P have the same origin: notably, the gf-kernel contributing to Ti [cf. f lA.4p - flA.5bp 
one appreciates from the relation 



as 
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y -In^ (y/y) 



ipn{x) + less singular terms in x. 



The crucial observation here is that the dispersion method yields an expression that contains 
the leading squared logarithm with a negative sign — in contrast to the result one finds with 



the factorization approach [cf. (13.130 ]. In this second case, it is more involved to show [44 
that the leading logarithm provides suppression in the relevant integration region and can 
therefore be associated with Sudakov effects. 



B. Spectral density in NNLO. /3o-proportional contributions 



The /3-dependent part of the partial spectral density, pn^\ can be obtained from the 
corresponding part of the whole amplitude T2, given by Eqs. flA.lQp . i.e., from Tp ® ipn- 
The calculation of the discontinuity of the latter expression is a rather technical task and 
is, therefore, relegated to the two Appendices \M and O It is important to realize that the 
structure of the spectral density, pn^^ , resembles the structure of the analogous term in the 
NNLO/3-amplitude that is proportional to T/^<^ipn, as one may appreciate by comparing the 
following two expressions: 



T2®tpn ho-Tp® tpniq^, Q^ pI, /Ur) = bo 



In { ^ ] Ti (g) 



2 



+ CpTo ® <; tP + Ln- 



(3.15a) 
(3.15b) 



rn 



bo ■ s^pI.pI) = boC^ [in (||) p«(^; s/pI) 



Rn\x;s/p 



(3.16a) 
(3.16b) 
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For convenience, we have made use of the abbreviation Ln [cf. (12. 7p ]. omitting arguments 
and separating for emphasis such terms from other functions by a dot. Recalhng the results 
of the previous subsection, one appreciates that the spectral density pn\x,s/ follows 
from the first term on the RHS of Eq. (I3.15al) . On the other hand, the new contribution 



Rn\x, s/ fip) in Eq. fl3.16bp derives from Eq. fl3.15bp and represents one of the main results 



of this investigation. 

To continue, consider the important partial case n = for which the corresponding 
expression for Pq^^"* becomes significantly simplified. Indeed, the leading logarithmic term 
contributing to R^^ and stemming from Ln^(?/) [cf. (]3.15bp ] cancels out because it is pro- 
portional to v{n). Actually, this is a general property of the leading logarithmic terms in 
all orders of the expansion that first reveals itself at the NLO level — see Eq. (I3.1al) . On the 
other hand, the subleading logarithmic term ~ In [s/np) survives, because it originates from 
the Vjj^^-element of the l^*-^-* kernel and from the T'^^^-NLO element of the hard-scattering 
amplitude that does not vanish for the ipQ harmonic [recall the discussion after Eq. (12. 8p 
Thus, the final result for rI^^ reads 



(x; s/i^l) = ® ^o) {x) + C2(x) + In ( ^ ) ■ C^{x) 



s 



(3.17) 



where the individual ingredients of this equation are the following 



Ci{x) 

C2{X) 



1 
'2 



■ipo = —6 (x ln(x) + X ln(a;)) 



du 



Ci(m)-Ci(x) 



+ (x x) 



3 -. 2/— \ 2"^ 

— xm (x) + vr — 
2 ^ ' 2 



u — X 

3 (xln(x) -|- X ln(x) — a;Li2(x)) + (x 



X) 



(3.18) 



(3.19) 



r^(2) = xx|30 [Li3(x) + Li3(x)] - 6 [ln(x)Li2(x) + ln(x)Li2(x)] - [in^(x) + In^(x)] 

-5 In' (I) + [ln(x) + ln(x)] (3 ln(x) ln(x) - 271^) - 72C(3) + ^tt' - 7} 

19 3 
+ — [xln(x) + xln(x)] + - [xln'(x) + xln'(x)] . (3.20) 

To derive Eq. (I3.17p . we have made use of Eqs. ( ]A.3[) . ( 1A.12P for the spectral density 
elements and also fl3.15bl) in conjunction with Eqs. flC.6p and fIC.lip . The coefficient Ci in 
front of In (s/pp) in Eq. (I3.17P accumulates those contributions responsible for the breaking 
of the conformal symmetry in the MS-scheme owing to V j45| and g j43| . 

Note that just this term leads to the breaking of the [x ^ x)-symmetry of the spectral 

density ^l^"* ( x, ) [recalling that s = Q'^x/x]. In Fig. [3l we compare the contributions 

to the spectral density from the NLO and NNLO/5 at one single scale Q' = /i| = /Xgy, the 
latter scale /Zgy = (2.4 GeV)' corresponding to the typical average momentum jsj measured 
by the CLEO Collaboration in the region [1.5 — 8] GeV'. This is done in terms of 
as{np)boCFR'o\x, x/ x) for the NNLO/3 spectral density (dashed green line) in comparison 
with p^o\x) calculated at the NLO (solid red line). 



13 




0.2 0.4 0.6 

FIG. 3: The dashed green hnc shows the expression as{^p) Pq^' = as(/i|) boCp R^\x , x / x) at the typical 
CLEO reference scale (labeled by the acronyms of Schmedding and Yakovlev [5[) 
(2.4 GeV)^, whereas the sohd red hne represents p'i^\x) in (|3.1ip . 



To render our presentation more transparent, we have relegated the discussion of the 
elements Tg'^\Vg^\T^^\ contributing to the amplitude T2 [see Eqs. ( I3.17p . flS.lSp ]. to the 



Appendices R] and iBl The expressions for the partial densities Rn\x , s / jji^) , appearing as 
convolutions of these elements with ifjm are supplied in Appendix [Cl 



IV. SIZE OF THE NNLO CONTRIBUTION TO THE FORM FACTORS 



In this section we discuss the effects of the NNLO/3 corrections to the transition form 
factors for the p and vr mesons, relying upon the BLM prescription and its modifications. 
Employing the results obtained in the previous section, we are going to simplify the expres- 
sions for the spectral density by adopting the so-called 'default' scale setting Pp = = Q^- 
Then, the expression for pn^\ given by Eq. fl3.16bp . reduces to 



pi^^)(x;.,g^g^) 



(4.1) 



(see Fig. [3]). Inserting this result in Eqs. (I2.24l) - (l2.25p and performing the integration over x, 
one finally finds for their sum [Eq. (I2.20p ] an expression that contains FpiQ"^) at the NNLO 
level. The final result for the whole form factor reads 



(4.2) 



It turns out that the NNLO/3 contribution, calculated here, is negative (dashed line in Fig. 3). 
Hence, taken together with the already known NLO contribution, which is also negative, the 
total effect of the radiative corrections at the considered level of the perturbative expansion 
is to decrease the magnitude of the form factor. Following the BLM procedure, the last term 
in Eq. (14.21) determines the "shift" of the scale in the argument of the running coupling from 
the value to the BLM-scale 42|, according to 



exp 



mi 



(4.3) 
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2 4 6 8 10 12 14 2 4 6 8 10 12 14 



FIG. 4: Left: The dashed (red) line in the middle represents Qblm(Q^)i whereas the solid (green) line 
corresponds to the result QIlm = exp(— 1.811) in [7|, and the upper dashed (blue) line marks the BLM 
scale for the V-part [cf. Eq. (|2.20p ] of the form factor describing the 7*p7r transition. Right: The upper 
dashed (red) line represents the BLM ratio i^blm, while the lower (green) line denotes 82 in Eq. (|4.7p . 



As a result, asiQ"^) — ^ cls{Qbi.m) > o.siQ'^) at <5blm(Q^) < hence, the form factor 

given by Eq. (14.21) assumes the BLM-improved form 

^^SrIQ') - F^iZiQ') = FoiQ') + as{QLM)F,{Q') . (4.4) 

The main contribution to -^lcsr provided by the asymptotic DA, ipQ. The associated BLM 

scale is QBLMlQmin) = ~ 1 GeV^ (see Fig.HJ left panel). This scale may be considered as 

the borderline for applying perturbative QCD, defining this way some minimal scale for the 

BLM scheme — denoted Q^j^. It is remarkable (though accidental) that, as mentioned above, 

this scale corresponds approximately to /igy = 5.76 GeV^ . Obviously, below this particular 

scale, the BLM prescription, expressed through Eq. (14.31) . would entail a renormalization 

scale that will be out of the region where perturbation theory can be safely applied. 

The scale Ql, determined above, still belongs to the perturbative regime and, therefore, 

the ratio 

SBLuifily) ^ [Flcllif^ly) " i^B^I^K/iiy)] / Flclnif^ly) ^ 0.11 (4.5) 

can serve as a crude measure for the relative weight of the NNLO contribution. In Fig. H] (left 
panel), we show Q^LuiQ'^) (dashed red line in the middle), following from Eq. (14. 3p . in 

comparison with a simpler linear dependence Qblm — exp(— 1.811) (solid green line) that 
results from the standard factorization formula of the perturbative approach employed in 0]. 
One observes that both results are rather close to each other in the important CLEO-data 
region. To be able to apply the BLM procedure below the minimal scale Qmin^ shall use 



a somewhat improved version of this procedure — termed BLM — introduced in (36[. Above 



the perturbation-theory borderline > Q^[^, this modified BLM procedure coincides on 
the RHS of Eq. (14.41) with the standard one. But for Q"^ < Qmin; the BLM scale is frozen at 
Ql and the expanded expression contains only a tail of the NNLO^ correction provided by 
the third term in the equation below 

F^iQ') = Fo(Q2) + a,(gjFi(g2) + aliQl)bo {Fp{Q') - F,{Q') HQ'/QD) . (4.6) 

Substituting this expression for F^^^ into Eq. (14. 5p . one obtains the quantity 5blm' shown 
in the right panel of Fig. |H The numerical value, estimated in Eq. (14. 5p . turns out to overes- 
timate the size of the NNLO contribution, as one observes from this figure. In mathematical 
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terms this becomes evident by glancing at the ratio 

6,{Q') ^ \al{Q')boF,{Q')\/F^'^i^{Q') (4.7) 

and recalhng that the magnitude of this contribution to the total form factor is a few times 
smaller than 5blm moderate values of characterizing the CLEO-data region. The 

size of the NNLO correction seems to be rather important and at the level of about 10% 
at low Q^, even by taking it into account only in the incomplete form of Eq. (14. 2p . On the 
other hand, at higher momenta, inspection of the 62 behavior in the right panel of Fig. H] 
reveals that the size of the corrections rapidly decreases to the level of 5% around the scale 



V. PREDICTIONS AND COMPARISON WITH EXPERIMENTAL DATA 

This section contains a discussion of the implications of our theoretical findings on the 
transition form factors F^*^'^ and F^''^'^ vis-a-vis the experimental data for the former. To 
understand the influence of the NNLO radiative correction on these hadronic observables, we 
have to analyze its relative weight with respect to the NLO contribution. Moreover, we have 
to discuss the interplay between perturbative corrections and nonperturbative ingredients, 
notably, the quark virtuality and the twist-four scale 5^. This discussion can be further 
substantiated by comparing the calculated photon-to-pion transition form factor with the 
available experimental data from measurements by the CLEO [9| and the CELLO lol 



Collaborations. From this comparison, we can extract valuable information as to what 
extent our calculation can describe the data in the whole measured region. From the 
theoretical point of view, we can use these data in order to estimate what is still missing 
on the theoretical side. We will focus below not on the exact phenomenological description 
of the mentioned data, but analyze instead the ramifications they imply on the theoretical 
approach and its various elements. 



A. Nonperturbative input 

We discuss first the nonperturbative ingredients of our analysis. The main one is the 
leading twist-two pion DA, , which can be derived with the help of various methods. 
These include — among others — QCD sum rules and lattice simulations. In addition, one has 



to model the twist-four component of the pion DA, see, e.g., the discussion in [19|, l21|, |22 
Lacking a complete derivation of the full pion DA from first principles of QCD, we are 
actually forced to reverse-engineer its structure from calculations of its first few moments 
(£^)^. To be more precise, one can calculate the moments of v^i^-* with standard QCD SR 



461 and also with those employing nonlocal condensates (NLCj fisl . 47, 48.j49l. 50|. On the 



other hand, one can use LCSR to analyze the high-precision CLEO data [9|]on F"'*"''^ and 
extract rather strict constraints on the first Gegenbauer coefficients a2 and [sl, 0, [l9| . More 
recently, two independent Collaborations have published results for the first coefficient 02 by 
measuring the first moment {^'^)tt of v^l^"* on the lattice jlsl, [H, [sH. The lattice calculation 
of the second coefficient ( or, equivalently, the fourth moment (^'^)7r) is still lacking, but a 
compatibility region between the CLEO data ^ and the a2-lattice constraints was worked 
out in jlsl to predict a rather narrow interval for the (^^)7r values. 
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FIG. 5: Left: Comparison of the CLEO-data constraints on '^JHQ'^) in the (a2, 04) plane at the scale 



. 1 

77r ( /^2\ 



fi'^Y in terms of error regions around the BMS best-fit point + @, Q, using the following designations: la 
(thick sohd green Hue); 2a (solid blue hne); 3a (dashed-dotted red line). Two recent lattice simulations, 



denoted by vertical dashed lines 



and solid ones 



29| are also shown together with predictions obtained 
from nonlocal QCD sum rules (slanted green rectangle) 18[ — the latter corresponding to the vacuum quark 
virtuality — 0.4 GcV^. Right: CLEO data in comparison with various theoretical models and lattice 
results listed in Table |I1 The dashed green la ellipse (together with the shifted best-fit point +) describes 
the effect of including the twist-four contribution to the pion DA via renormalons [22 1 . 



We provide below a short overview of our present knowledge of (f n from different sources, 
omitting specific details for which we refer the interested reader to the original literature. 
The models shown below in the figures are summarized in Table [H Using QCD sum rules 
with NLC, we derived a "bunch" of admissible pion DAs [3], taking into account in the 
expansion fl2.10p only the first two terms with 02 and 04 (details can be found in [s^). Pion 
DA models, like the Chernyak-Zhitnitsky (CZ) model [i^, or the Braun-Filyanov (BF) 



38l | one, also used these two harmonics for modeling the pion DA. Note, however, that in 
the NLC approach [3] this is not the result of an arbitrary truncation of the Gegenbauer 
expansion after n = 4, but follows from the fact that all calculated higher-order coefficients 
up to n = 10 turn out to be compatible with zero. Hence, from a pragmatic point of view, in 
order to capture the main characteristics of the pion DA, it is sufficient to restrict the analysis 
of the experimental data on the photon-to-pion transition to two-parameter models. Such a 
data analysis of the CLEO measurements was first performed by Schmed ding and Yakovlev 
(SY) 1^ within the framework of LCSR. In a subsequent series of papers [6|,[U 2^, 25, 26, 35 
(nicknamed BMS), this type of data processing was further pursued with results confirming 
the previous SY findings while also improving the theoretical accuracy. These results are 
displayed in Fig. [5] with details being provided in Table [H In this graphics (left panel), the 
slanted shaded (green) rectangle represents the BMS DA "bunch" from the NLC approach 
together with the middle point X, which corresponds to the BMS model with the coefficients 
a^^^inly) = +0.14, a^^^d^ly) = -0.09. The CLEO-data constraints are shown in the form 
of error ellipses: la (thick solid green line), 2cr (solid blue line), and 3a (dashed-dotted 
red line). The range of values of a2, determined recently on the lattice by two independent 
Collaborations, is denoted by vertical dashed lines [2^ and solid ones 2§]. All constraints 



and predictions shown have been evolved to the scale /igy using the NLO ERBL evolution 
equation. 

In the right panel of Fig. [5] we show the effect of the twist-four contribution to the pion 



DA, taken into account via the renormalon approach of [60| , in comparison with the standard 
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TABLE I: Estimates of the Gegenbauer coefficients 02 and 04 at the normahzation scale //gy 
obtained by different methods and for several pion DA models. The designations correspond to 
those used in Fig. [5l Also included are the theoretical constraints derived from QCD SRs with 
NLC and estimates derived from an analysis 0] of the CLEO data 0] using LCSRs. The lattice 
measurements of [29|, [28[, and 5l|] are also shown. [Note that the uncertainties of the coefficients 
02 and 04 are correlated. Here, the rectangular limits of the fiducial ellipse [g, 19| are shown.] 



DA methods/models 


Symbols 


«2(AtsY) 




Lattice 








UKQCD/RBC [29] 


1 1 1 


0.215 ±0.07 


— 


QCDSF/UKQCD [28] 


[ 1 


0.19 ±0.11 


— 


[5]J 




0.233 ±0.145 


— 


Data analysis 








LCSR for FT", CLEO data 








SY best fit [5] 


• 


0.19 ±0.04 ±0.09 


-0.14 ± 0.03 =F 0.09 


EMS best fit [6, 19} 


+ 


0.22 


-0.22 


BMS \9,19}{la) 




[0.11, 0.328] 


[-0.03,-0.41] 


Guo&Liu \27}{la) 




[0.06, 0.14] 


[-0.02,-0.13] 


BMS best fit [22] 








(twist-four via renormalons) 


+ 


0.31 


-0.25 


Agaev best fit [52] 


A 


0.23 


-0.05 


LCSR for F'", JLab data. 








BK model [53 




0.17 ± 0.11 




NLC QCD SR 








QCD SRs for (^^)^ [18, 35] 




[0.1, 0.185] 


[-0.03,-0.14] 


BMS model [18, 35] 


n 


0.14 


-0.09 


Instanton models 








ADT 1541 


▲ 


0.034 


—0.027 


PPRWG [55] 




0.03 


0.005 


PR [56] 




0.06 


-0.01 


NK [57] " 




0.119 


0.014 


Models 








As (Asymptotic) 


♦ 








BZ [58] 


▲ 


0.08 


-0.01 


CZ [46] 


■ 


0.40 





BF [38] 


T 


0.31 


0.15 


AdS/QCD [59] 








(a2, 04) projection 


• 


0.1 


0.035 



"This instanton model improves the works in 
CLEO data. 
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56l ] and turns out to be inside the 2a error ellipse of the 
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one, which is based on the asymptotic DAs of twist four namely, 



V>^^\x,f,') = ^-^6\f,')x\l-xr. (5.1) 



Here 5 (IGeV ) = (0.19 ± 0.02) GeV . This estimate has been obtained in [6|] under the 
assumption that = {q {ig a^^G^"") q)/{2{qq)) = 0.4 ± 0.05 GeV^ (see Appendix A in Q). 
A full-fledged analysis of the renormalon-model corrections has been given in [l^l — see also 



2l| . The net effect is to shift the whole la error ellipse (broken contour in Fig. [5]) further 
away from the asymptotic pion DA and hence to higher values of a2 (observe the shift of the 
best-fit point to its new position denoted by +). The original la error ellipse (solid contour) 
is also shown in this figure and one appreciates that the variation of the twist-four term can 
have substantial influence on the transition form factor. 

Let us conclude this subsection with the following remarks: (i) The most striking message 
of Fig. \5\ (to be read in conjunction with Table [I]) is that the CLEO la error ellipse, the 
two recent lattice calculations, and the fiducial region of pion DAs extracted from nonlocal 
QCD sum rules all have a common region of validity. Moreover, the BMS model DA K is 
entirely within this area, (ii) In contrast, the asymptotic pion DA ♦ is outside the 3a ellipse, 
while the CZ pion DA ■, as well as the BF model ▼, are outside the 4a error ellipse of the 



CLEO data, (iii) All models inside the "rhombus" , determined in [58| , are also more or less 
outside the la error ellipse. This fact gains more weight in view of the reduced interval of 02 
computed on the lattice j29| (in the form of {^^)n) which seems to exclude all these models 



as well, (iv) Analogous considerations for the moment (^^)7r have been discussed in 26 



B. Phenomenological input of the light-cone sum rules 

The next task concerns the model of the resonances entering the dispersion integral for 
the transition form factor. As we mentioned before, we are going to refine the simple 6- 
function ansatz by taking into account a finite width in terms of the Breit-Wigner model. 
To this end, we replace in 



pP'^-(s, = 0(^so - s)p\s, Q') + eis - so)Nr p{Q\ s) (5.2) 



by 



V 2vr ^t;;^ (m^. - s) + mlVl 

with ~ /p/3,and p^'^'^ ~ SF^''^'^. Then, we trade the simple 5-function resonance 
model for the Breit-Wigner ansatz [13], given by 

5{s - ml) ^ ^y{s) ^ \ ^ "^f^ ^ , (5.4) 



also taking into account the difference in the masses and widths of the p and the io vector 
mesons, rup = 0.7693 GeV, = 0.7826 GeV, and = 0.1502 GeV, F^ = 0.00844 GeV, 
respectively. To continue, in order to obtain F'^ , we appeal to the duality between the 
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hadronic part of the spectral density and its perturbative counterpart that we express via 



'""'■('• e'^.^ivj^^i.. (5.5) 



J s + q 



Ami. 



The main effect of using Eq. (15.41) on the LHS of Eq. (15.51) is a shght suppression of the 
integral relative to the outcome of the simple (5-function ansatz. 

After the Borel transformation of Eq. (15.51) . one arrives at an expression for F'^ in 
terms of V [cf. f lOS]) ]. notably, 



m = k~'^V{Q\M') , (5.6a) 



So 



Ap(s) + A, 



4ml 



where k ~ 0.932 weakly depends on the Borel parameter in a region containing the 
standard scale ~ 0.7 GeV^. Therefore, the Breit-Wigner ansatz for the resonances 
(abbreviated in what follows by BW) supplies a factor > 1, entailing an increase of the 
value of the form factor as compared to the simple (5-function ansatz. A similar enhancement 
effect in analogy to Eq. (12.201) applies to the total form factor -^lcsr' leading to 

F^clliQ') = Nt {I^^^W''^') + ^^(^')} ' (5.7a) 



so 



k,= / ^(£)+A(£) , (5.7b) 



where ki ^ 0.984. Crudely speaking, the net enhancement amounts to 2-4% in that mo- 
mentum region, where the resonance part prevails over the pointlike one. 



C. NNLO effects on Q^F'^*'^'^ — integral characteristics 

Having set up the framework for calculating the pion-photon transition form factor, in- 
cluding also its nonperturbative ingredients, we now turn our attention to the specific effects 
due to the NNLO/3 radiative corrections according to Eq. (14. 2p . 

Being interested mainly in the magnitude of the form factor, it is actually sufficient to 
include only the ipo term. An analysis including higher harmonics will be presented in a 
future publication. 

The final results for the photon-to-pion transition form factor, including the NNLO^ cor- 
rections, within the improved LCSR approach are displayed in Figs. and [71 The presented 
results have been obtained along the hues followed in our previous dedicated and detailed 



works in [6|, Il9|, l22[ and will not be repeated here. Recalling the features of the pion DAs 
discussed in Subsec. IV At we reduce our discussion to only three models: As (lower black 
line), CZ (upper red line), and BMS (green line in the middle). The solid lines represent the 
scaled form factor with the NNLO/3 and the BW-ansatz corrections included, whereas the 
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FIG. 6: Theoretical predictions from LCSR for the form factor Q^F^ '''^ with the NNLO^ corrections 
(together with the BW resonance model) included (solid lines) and without them (dashed lines) . All predic- 
tions shown are evaluated with the twist-four parameter value a = 0.19 GeV^ @, [l^. They correspond 
to selected pion DAs. These are: CZ model — upper (red) line 46|, BMS-model — middle (green) line 18 1, 
and As DA — lower (black) line. For comparison, the corresponding predictions for each pion DA model 
without these corrections are displayed as dashed lines. The experimental data shown are from the CELLO 
(diamonds, [l^l) and the CLEO (triangles, Q) Collaborations. 



dashed lines (with the corresponding color for each model) show the result without these cor- 
rections. These theoretical predictions are displayed in the background of the experimental 
data from the CLEO (triangles) and the CELLO (diamonds) [13] Collaborations. 

The main lesson from these figures is that the effect on Q'^F'^ '^'^ , induced by the radiative 
corrections, amounts to about (—10 -. — 5)% (see Sec. [IV]), whereas that caused by the 
use of the more realistic BW-resonance ansatz provides a small growth of approximately 
(+4-=-+2)%. Combining both effects, results into a net reduction of the magnitude to within 
7% at small — 2 GeV^. In the case of the As DA, the size of this reduction rapidly drops 
to 2.5% for Q"^ > /igy , as one may appreciate by comparing the solid line (NNLO/3) with 
the dashed line (NLO) at the bottom of Fig. [61 These results have been obtained with the 
spectral density p^Q^\x] s, Q^, Q"^) = C-pRq'\x, x/x), where the further evaluation goes along 
Eqs. (I3.18p - (l3.20p for the -j/^o^harmonic. On the other hand, the transition form factor is 
mainly determined via the inverse moment jlH, which belongs to the integral characteristics 
of the pion DA, i.e., 

{x-^)^= [ dx^^ = 3(l + a2 + a4 + ...) , (5.8) 
Jo ^ 

that can be easily understood from the expression for the H part in Eq. (12.251) setting 
3> So, — »■ 1. This means that for pion DAs, like the BMS one, which have the 
particular property 02 ~ — 04 (see the anti-diagonal in Fig. [5]), the corrections associated 
with the higher harmonics ip2 and ip4 mutually cancel, leaving only the correction due to ipQ. 
The same conclusion can be drawn for all pion DAs belonging to the BMS "bunch" — shown 
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as a green strip in Fig. [3 — (consult Sec. IV Al for further explanations).^ This cancelation 
effect is absent for the CZ DA, as one observes by comparing the solid (red) line (NNLO/3) 
with the dashed (red) line (NLO) on the top of both these figures. The reason is that the 
CZ DA has only a large 02 coefficient, while all other coefficients are zero. Hence, we cannot 
estimate the size of the NNLO correction using only the ipQ part of the NNLO contribution, 
but we have to include the corrections for higher harmonics as well. 
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FIG. 7: The same designations as in the previous figure apply. The shaded (green) strip denotes the 
form-faetor predictions derived with the help of the BMS "bunch" 18|. The dashed blue line represents a 
dipole-form interpolation formula [ol] of the CLEO data. 



We observe from Figs. [6] and [7] that the improved theoretical predictions which contain 
the discussed effects fail to describe both sets of the experimental data at low Q^, say, below 
4 GeV^. Actually, the data can be grouped into three regions: 

(i) The low region is mainly covered by the CELLO data and extends to momenta up to 
approximately 1 GeV^. This region, where hadronization is immanent, is virtually inacces- 
sible to the methods used in our analysis — let us call it therefore the unknowable regime. 

(ii) There is some intermediate momentum region between 1 and 3 GeV^, where the existing 
data are underestimated by our theoretical predictions. In contrast to the previous case, 
here we can figure out what the origin of this drawback may be. Hence, let us call this inter- 
mediate domain the unknown regime, because the missing contributions can, in principle, 
be estimated. 

(iii) Finally, above about 3 GeV^, the agreement between the CLEO data and the shaded 
BMS strip is fairly good, as one also realizes by comparing these data points with the dipole 
fit (dashed blue line) used by the CLEO Collaboration. In fact, the dipole curve and the 
prediction due to the BMS model (solid curve within the BMS strip) almost coincide in this 

region. Therefore, it is safe to claim that this high-momentum domain is well-reproduced 
by our techniques — hence the term known regime.^ 



^ Note that the width of the BMS strip is somewhat narrower compared to our previous results in 35 1 
because here we use a smaller range for the 5^ uncertainties in the twist-four eontrib iitio n. 

^ However, one may face the challenge posed by the recently released BaBar data 6l| exactly in this 
momentum-transfer region — see Note Added below. 
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Let us now discuss the reasons for the discrepancy in the unknown regime. 



One reason for the observed discrepancy may be traced to the fact that we used for 
the evolution of the twist-four contribution only the one-loop anomalous dimension. 
Considering the evolution effect at the two-loop level, could potentially reduce its size, 
thus rendering the reduction of the transition form factor less severe. 

The uncalculated remnant of the NNLO contribution could eventually enhance its 
total magnitude — should this part appear with the opposite sign with respect to the 
calculated /^-function part. 

Still another source of uncertainty lies with the value of 5^ which controls the size 
of the twist-four contribution. A smaller value of this parameter would entail less 
suppression of the form factor, given that this nonperturbative contribution has a 
negative sign, just like the calculated NNLO/5 contribution. Here there is a subtlety. 
If one assumes a smaller value of 5^, then this would automatically mean that also 
\g should assume a smaller value because 6'^ ~ -^5/2, i.e., these two nonperturbative 
parameters work synergistically. But a decrease of the latter parameter would yield to 
an enhancement of the leading-twist contribution, meaning that the whole BMS strip 



(made up on the basis of the BMS "bunch" [18|, |35|) would move somewhat upward 
as a whole. As a result, one would obtain an increase of the transition form factor 
exactly in that regime between 2 and 4 GeV^, providing a better agreement with the 
CLEO and CELLO data. 



D. NNLO effects on Q'^F'f* {Q"^) — differential characteristics 

The scaled form factor Q^F'^*p'^{Q'^), obtained in the framework of LCSRs, see Eq. fl5.6l) . 
is determined by the hadronic part V, defined in Eq. f l2.24l) . At moderate around the 
scale /igY, V is mainly formed by the y^^-dependent leading-twist contribution and to a 
lesser extent by the twist-four one. The combined effect of the BW-ansatz and the NNLO/3 
corrections is rather important and amounts to approx. —10% in the region limited from 
above by the scale /Xgy. For still higher momentum values, and up to ~ 15 GeV^, it 
reaches the level of +9% — see left panel of Fig. 8. More specifically, the NNLO/3 corrections 
reduce the value of the form factor by an amount of approx. 20%, starting at = 2 GeV^, 
and become significantly smaller at the end of this region (see dashed line in Fig. 8, left 
panel), whereas the use of the BW-ansatz leads to an overall increase of 8%. The ratio 
of the combined effect of the NNLO/5 and the BW-ansatz contributions, relative to the 
total transition form factor, is illustrated in Fig. [8] in terms of a solid red line, whereas the 
normalized contribution of the NNLO/3 term alone is depicted in the same figure by the 
dashed blue line. 

The form factor Q"^ F'^* (Q^) is presented in Fig. 8 (right panel) using different pion DAs. 
We observe from this figure that this process can actually be used to discriminate among 



different pion DAs of twist two, as first pointed out by Khodjamirian in [IJ] because it is 
sensitive to the particular shape of the pion DA. Especially the momentum region between 
2 GeV^ and 8 GeV^ seems to be particularly convenient for this task because (i) the obtained 
predictions are clearly distinguishable and (ii) because it corresponds to the range probed 
already by the CLEO experiment for the pion-photon transition [9|. The measurement of 
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both transition form factors in this momentum region would provide definitive clues for the 
underlying pion DA. 
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FIG. 8: Left: The upper solid curve shows the combined contribution of the NNLO/3 and the BW-ansatz. 
The dashed hue below it represents the normalized contribution of the NNLO/3 only. Right: Predictions for 
Q4^7Pir^Q2-j jtqj, different pion DAs: the top solid (red) line denotes the CZ DA, the middle (green) strip 
represents the BMS "bunch" — including the BMS model DA as a dashed line — , and the lower solid (black) 
line gives the result for the As DA up to 15 GeV^. The other lines shown inside the BMS strip are explained 
in the text. 



Let us now make some detailed remarks on the structure of the BMS strip, obtained from 
the "bunch" of pion DAs following from NLC QCD sum rules. The first observation from 
Fig. 8 (right panel) is that the BMS strip (which includes the BMS model: long-dashed 
line) [l8| crosses around 11 GeV^ the prediction obtained with the As DA. We argue that 
this effect can be traced back to the differential pion characteristics of the pion DA, viz., 

V3^(0) = -j-(Pn{x)\x=o- The usefulness of the slope of the pion DA was first discussed in H 



dx 

in connection with the pion's electromagnetic form factor. In our case, it is instructive to 
recall the definition of V in Eq. (I2.24p . from which one appreciates that the lower limit of 
the dispersion integral xq tends to the upper one, 1, when ^ sq. Therefore, near the 
upper integration limit, this integral is determined by the endpoint behavior of p{x)\xn^i, 
which in LO is proportional to (/'^(O) [cf. Eq. (I2.9al) ]. entailing for the form factor a behavior 

like ~ V^^(O) 2- This becomes visible in the vicinity of the upper limit 1, say, for 

(Q^ + So) 

a value 0.1 = sq/{sq + Ql), which for the BMS model corresponds to Ql ~ 14 GeV^, while 
for the upper part of the strip the value of Ql is much larger. Thus, it becomes pos sible 



m 



to investigate the endpoint behavior of different pion DAs — mentioned also in [IJ, |38 
terms of this transition form factor by probing the slope of the pion DA at the origin. 

To make these statements more quantitative, let us express the slope of the pion DA in 
terms of the Gegenbauer coefficients to obtain 

(^;(0;/i2) = 6 (1+ {m+l){m + 2)/2aUfi^)] . (5.9) 

\ m=2A,... J 
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Then, we get the following results: 



17.4 [CZ] 

^'MQl)= { 3.7 [BMS] (5.10) 
5.8 [As] . 

These simple estimates are qualitatively responsible for the associated form-factor values 
at (slopes and altitudes) in Fig. [8] (right panel). It turns out that for the CZ DA the 
value of the form factor is approximately 3 to 4 times larger than the one for the As DA, 
while for the BMS DA it is smaller or quite close to it [cf. Eq. (15.101) ]. Note that this 
feature appears to be opposite to the Q'^F^*^'^{Q'^) form factor that mainly depends on 



the inverse moment 22| {x~^).,^, the latter being an integral pion characteristic [18|, [l9 
Indeed, one sees from Fig. [7] that the results for Q'^F^*^'^ {Q"^) (green-shaded strip) are 
always larger than the predictions obtained from the asymptotic DA. Hence, Q^F'^'^'^IQ'^) 
can provide complementary information about the pion DA and help discriminate among 
various proposed pion DA models. 

Going beyond the leading order of p, a new behavior near the endpoints emerges due to 
the hard-gluon exchange leading to a 1 / Q^-behavior of the form factor in the far asymptotic 
domain. This is, because in NLO, the expression p^^^ in Eqs. (13.111) and (13.121) contains a 
term proportional to —ln^{x/x) that accumulates this effect. 

A second issue related to the BMS strip in the right panel of Fig. [HI which deserves to 
be considered as well, concerns the fine structure of its envelopes. Indeed, close inspection 
reveals that the upper envelope has a dip at ~ 9 GeV^, whereas the analogous irregularity 
of the lower envelope is much less pronounced and appears at a scale close to 10.5 GeV^. The 
origin of these irregularities might be related to the crossing of the predictions obtained from 
different DAs inside the BMS "bunch" that have different initial values of the Gegenbauer 
coefficients 02 and 04 (exemplified by the intersecting lines inside the BMS strip). When it 
happens that such a crossing point lies just near the upper or lower boundary of the BMS 
strip, the corresponding envelope is "bent" inwards. Physically, this means that the range 
of the calculated uncertainties is somewhat smaller there relative to the regions where no 
crossing points appear close to the boundaries. On the other hand, we cannot exclude that 
these irregularities are a spurious effect created by the algorithm we used to obtain our 
predictions. In any case, they have no influence on our analysis. 



VI. CONCLUSIONS 



Results have been presented for the photon-to-pion and the p-meson-to-pion transition 
form factors using light-cone sum rules and including those NNLO contributions which are 
proportional to the /5 function. Let us summarize our main flndings. 

1. The spectral density of the LCSR at the NLO level was systematically constructed for 
any Gegenbauer harmonic of order n and was presented in compact form in Eqs. (13. ip 
and (Km . 

(2) 

2. Using the NNLO hard-scattering amplitude Tg , proportional to the (3 function and 
calculated before in 0], we derived all necessary ingredients of the NNLO spectral 
density pn and analyzed Po for the harmonic via Eqs. (I3.18I) - (I3.20I) . These 
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quantities enter the dispersion integral which determines the structure of the form 
factors '"^ and F'^ ^'^ within the hght-cone sum-rule approach. 

3. Predictions for the form factor F^*"^^ were presented which include the NLO and 
the /3-part of the NNLO radiative correction, an improved resonance contribution, 
and also estimates of the uncertainties stemming from the twist-four contributions. 
The first contribution has a negative sign and amounts to a reduction of up to 10%, 
whereas the Breit-Wigner ansatz provides a small enhancement below 5%, and the 
twist-four contributions can yield to a small enhancement by adopting a smaller value 
of the parameter 5^ which controls their size. This, however, would entail a parallel 
enhancement of the leading-order contribution to F^*^'^ via the parameter A^, because 
these two parameters are correlated (5^ ~ Ag/2). Comparison with the CLEO and the 
CELLO data shows that at 2 GeV^, and below, the calculated form factor falls slightly 
short compared to the data, leaving room for additional soft and higher perturbative 
contributions that are unknown at present. 

4. We have also given predictions for the F^''^'^ form factor using different pion distri- 
bution amplitudes including the asymptotic one, the whole set of pion DAs derived 
from nonlocal QCD sum rules, and the CZ model. In this case, the combined effect 
of the negative NNLO^ radiative corrections and the use of the Breit-Wigner ansatz 
lead to a net reduction of the form-factor magnitude varying between —10%, below 
the scale /igy, and growing to the level of +9% at ~ 15 GeV^, whereas the twist- 
four contributions amount to —11% even up to the scale 15 GeV^. Once there will 
be experimental data for this reaction, it will provide an additional and useful tool to 
discriminate among various pion DAs — especially between those of the CZ type that 
receive strong endpoint enhancement on the one hand and such with their endpoints 
being suppressed — like the BMS one — on the other. Moreover, because this form factor 
is sensitive to the differential characteristics of the pion DA — expressed via '/'^(O) — one 
can use the endpoint behavior of the underlying pion DAs as an additional adjudicator 
in selecting the optimal pion DA. 

To conclude, we have extended the analysis of the pion-to-photon transition process to 
the NNLO level in a systematic though partial way. Nevertheless, our analysis provides the 
possibility to estimate the size of the associated form factors in that region of momenta 
which is accessible to present measurements. It will be interesting to pursue and complete 
this sort of calculation in the future by including the whole NNLO contribution. 
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Note Added 



After completion of this work we became aware (thanks to Dieter Miiller and Maxim 
Polyakov) of new prehminary data of the BaBar Collaboration^ on two-photon-induced 
processes in the momentum range 4 < < 40 GeV^ which show above 10 GeV^ "a power- 
law growth behavior that contradicts most models for the pion DA" . Meanwhile, these data 



have been officially released [61 . 

Taking the high-Q^ data points of BaBar at face value, one might even come to the 
conclusion that they are incompatible with QCD per se because the indicated power-law 
enhancement at large momentum values cannot be explained by any known QCD effect, like 
higher-order radiative corrections or contributions due to higher-twist effects. Our presented 
analysis serves to prove that the inclusion of the main part of the NNLO radiative corrections 
provides suppression — not enhancement — at the 10% level, while the twist-four contribution 
gives a small enhancement of a few percent in the CLEO region. The expected size of the 
uncalculated NNLO remainder — even if it should have a positive sign — is not expected to 
exceed the few-percent level (< 10%). An enhancement of the twist-four contribution due 
to two-loop evolution is possible but it is expected to be of the order of a few percent as 
well. The size of all these QCD corrections is far less than indicated by the high-Q^ BaBar 
data, so that the observed enhancement at high cannot be explained by higher-order 
perturbative QCD and power corrections. 

Note that endpoint-enhanced pion DAs, like the CZ model, are also in conflict with these 
data, despite opposite claims in [61[, because the corresponding Q'^F'^''^'"{Q'^) prediction (see 
Fig. inj scales with (analogously to all other pion DA models in the convolution scheme) 
above approximately 15 GeV^, in sharp contrast with the signiflcant growth of the BaBar 
data in this momentum range. In this context, we flnd it remarkable that the BaBar data 
points in the momentum range already probed by the CLEO Collaboration are compatible 
with the CLEO dipole flt and the asymptotic QCD prediction v^/vr, being also within the 
BMS strip. Even more signiflcantly, two more data points — outliers — at about 14 GeV^ 
and 27 GeV^ turn out to be just on the upper boundary of the BMS strip and, hence, in 
compliance with the QCD expectations, (see Fig. [9]). Hence, one may divide the BaBar data 
into two branches: one containing the data points in the CLEO region plus the two outliers, 
the other consisting of the remaining 10 high-Q^ data points. The flrst branch supports 
the QCD predictions with NNLO radiative corrections and twist-four contributions. The 
other branch is in clear conflict with the convolution scheme of QCD. In view of this data 
structure, odds are that the BaBar data may bear some intrinsic inconsistency. 

We end this discussion with the following key observations: 

(i) The main NNLO radiative corrections and the twist-four contributions do not provide 
enhancement to Q'^F"'*"''^{Q'^) in the range of momentum transfer 10-40 GeV^, exclusively 
covered at present by the BaBar experiment. Hence, the observed behavior of Q'^F^'^'^iQ'^) 
growing with Q"^ in this momentum region with an almost constant slope cannot be ex- 
plained within the convolution scheme of QCD. 

(ii) Within the QCD convolution scheme, all pion DA models, which have a convergent 
projection onto the eigenfunctions (Gegenbauer polynomials) of the meson evolution equa- 
tion and hence vanish at the endpoints and 1, are conflicting with the BaBar data for 



^ Talk presented by Selina Li on behalf of the BaBar Collaboration at Photon 2009, Hamburg, Germany, 
ll-14th May 2009. 
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TABLE II: Deviation in terms of = x^/ndf (ndf = number of degrees of freedom) of Q'^F^ ^^{Q"^) 
predictions for the asymptotic (Asy), the BMS DA, and the CZ one. For a direct comparison with 
the BaBar analysis, we employ the NLO approximation of perturbative QCD, including also twist- 
four contributions. The first column shows the results for the combined sets of the CLEO and 



the BaBar [6l[ data. The second column refers only to the BaBar data, while the third column 
takes into account only the last 10 high-Q^ BaBar data, starting with the data point at 10.48 GeV^. 



Pion DA models 


CLEO and BaBar 


BaBar (all data) 


BaBar (10 data > 10 GeV^) 


Asy [31, 32j 


11.5 


19.2 


19.8 


BMS [18j 


4.4 


7.8 


11.9 


CZ [46J 


20.9 


36.0 


6.0 



g2^7*7^(g2) between 10 and 40 GeV^. 

(iii) Staying within this approach, the best agreement to the combined sets of the CLEO 
and the BaBar data is still provided by the BMS-type pion DAs, as one sees from Table [TTl 
first column. Considering only the BaBar data, the deviation of the CZ DA becomes even 
larger — second column in the same table. The CZ DA is favored only when one includes in 
the fit the last 10 high-Q^ data points, starting at 10.48 GeV^ (third column in Table HTl) 
But such a treatment of the existing data (CLEO and the 7 lower-Q^ BaBar data points) 
looks biased or at least unjustified. From this discussion it becomes clear that the conclusion 
drawn by the BaBar Collaboration in that the CZ DA is in agreement with their data 
(within the convolution scheme) is unfounded. 
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FIG. 9: Predictions for (52^7*7^(^2^ derived from three different pion DA models: Asymptotic, BMS, and 
CZ. The same designations as in Fig. [71 are used. The momentum-transfer range is extended to 40 GeV^ in 
order to include the new BaBar data [611 (shown as thick bullets with error bars). The displayed theoretical 
results include the NNLO/3 radiative corrections and the twist-four contributions. The horizontal dashed 
line represents the asymptotic QCD prediction v^/tt- 



The bottom line: The anomalous high-Q^-behavior of the pion-photon transition form 
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factor, found by BaBar, could potentially be of great importance, if confirmed by indepen- 
dent measurements, e.g., by the BELLE Collaboration, and would demand a new framework 
of analysis within QCD. Some scenarios to explain the BaBar effect have already been pro- 



posed [6^, l63|, l64 



APPENDIX A: DETERMINING THE SPECTRAL DENSITY FROM THE DIS- 
CONTINUITY OF THE HARD-SCATTERING AMPLITUDE 

Let us first define the task ahead in a bottom-to-top approach, starting with the LO 
case and finishing with the NNLO one. First, we analyze the structure of the amplitudes 
themselves, then, we construct the discontinuity for their elements in order to pair each 
Tk with the spectral density, p^''^ , associated with this order of the perturbative expansion. 
At the LO level we encounter just one single contribution to the discontinuity, i.e., that 
one arising from Tq ® / (where f{x) is some appropriate case- dependent test function). In 
NLO, Ti ® /, a new element appears: Tq ® Ln(a;) ® f{x). Finally, at the NNLO level of 
the expansion, T2, we are faced with the discontinuity of still one more quantity, notably, 
Tq ® Ln^(a;) ® /(x). [Recall that Ln(x) denotes the logarithm of the photon momenta over 
the factorization scale — cf. fl2.7l) .] 

Before we continue, we mention parenthetically in this context that the contributions 
to p^''\ just mentioned, can be also obtained from a generating function R(s;e), which we 
display below: 



R{e, x;s,Q ) = ' ' = — — ^ x ^ (x - xq 



s + g2 2(s + g2) V /iF ; ' ' tt 

+{x x) (A.l) 
P^'\Q\xo) = (^ly i?(e,x;.,g2)[^^. (A.2) 

The total sum of these contributions, R{as,x] s^Q"^), (following from Eq. ( 1A.1I) ) enters the 
spectral density in the form of a convolution with R® Therefore, the resummed 
expression does not contribute to the ipo part by virtue of the current conservation f (0) = 0. 
1. From the discontinuity of Tq in LO 

To(g^ -S] x) = ^ ( TTT- + ^ X 



we have jl, [l3l 

p(^\Q\s) = -Im 



2 \Q'^x — sx 



ToiQ , -s) ®(p^= / + (x ^ x) (p^{x)dx 

^{Q + s) Jo \X-Xo 



2(Q2 + s 



6{x — Xo) + (x — > x) )ipT,{x)dx 

(A.3) 



Xo = s/{Q' + s) 



2. In NLO we obtain 0, Q 

Tl(g^gVF;a:)=CF|^«(g^g2;x) + ln(^^^ [Tq ® y+] (g^ x) | , (A.4) 
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T^^\Q\q^-x)=To{Q\q^;y)®[T^^\y,x) + m{Q\q^-y)V{y,x)+\ , (A.5a) 
r^^\x,y)=[-W' + g] (x,y)+-35(x-y), (A.5b) 

where the notations of Ref. [3] and the following abbreviation have been used: 

Here Q'^ = — [(gi — q'2)/2]^ = {Q^ + ?^)/2, and ii\ and yUp denote, respectively, the scale of 
the renormalization of the theory and the factorization scale of the process. The kernels 
and V'^ are diagonal with respect to V'n and are defined in Appendix [HI whereas the kernel 
g reads 



LN(Q2,g2.^)=ln 



(A.5c) 



9ix,y)=-2 



0{y - x) 
y-x 



In ( 1 - - 



(A.6) 



and is not diagonal with respect to the Gegenbauer-polynomials. This kernel is responsible 
for the apparent breaking of conformal symmetry in the MS-scheme 0] • Substituting flA.Sap 
and flA.ScI) into Eq. flA.4p . we arrive at 

Tl(g^ q'; fil; x) =UQ\ q'; y) ® {C^T^'\y. x) + Ln(y)l^W(i/, x)+} . (A.7) 

One observes that in this expression a simpler logarithm Ln(?/) = In [{Q'^y + q^y) / ap- 
pears in comparison to LN in Eq. (]A.5cp . From the former expression and Eq. (lA.5bp . we 
get the following result for the convolution 



Ti ® ^„ = C^To{Q\ g'; y)®[ - AM ■ 3 [l + v\n)\ + [g+ ® A] {y) 

+ ^Jnijj) ■ 2v{n) Ln{y) 



where the second term in Eq. (lA.SP reads 
[g+0i;n]iy)= ^-InM^ 

L 3 \y 



1 r 



) - i^niy) 



In 1 - 



x-y 



x-y 



In ( 1 - - j dx+{y ^y) 



dx+{y y) 

^ Gni My) ■ 



(A.9) 



1=0,2.. 



Here we present a few partial values of the Gni elements: 



Goo 
G21 

Gqi 



3 35 
2'~12 
3 161 

2' ~ri 



203 



45 

3 49 781 



180' 40' 300' 



29531 1 
5040 J ■ 



(A.IO) 
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The RHS of Eq. (1A.9I) can be further evaluated using the following relation 



ry rj.1 _ yl 

lo y-x 



ln(l 



dx = y 



+ + . (A.ii) 



The discontinuity of Ti ® ipn in Eq- (lA.lOl) consists of one part, determined by the disconti- 
nuity of To and following directly from Eq. flA.Sp . and a second nontrivial part entailed by 
the logarithm, as one can also verify from the analysis in j^. Then one has 



Im 



In 



iQ' + s)y 

i4 



1 



dy 



2(S + y-y^ 
1 



(ln(y - yo) - In(yo)) + In ( — 



yo) 



2(g2 + s) 
+iyo yo 



dy 



d{y-yo)\n { — 
Pp 



djyo > y) 
y-yo 



Hyo) 



(A.12) 



Substituting Eq. flA.Op into flA.lOl) by taking into account Eqs. f lA.Sp and flA.121) . and 

collecting all terms, one finally arrives for the spectral density Pn\y) at the final results 
given by Eqs. fl3.lap - fl3.lcl) . Note that the last term in Eq. flA.12p generates the term f]3.1bp . 
To obtain the "algebraic" form fl3.12p . one should insert the relation (13.50 into the LHS 
of Eq. (13. 7p . The first partial hni coefficients of the expansion of the final result over the 
Gegenbauer harmonics are determined to be 



^00 

h2i 
hi 



3 

2' 



3 


25] 




2' 


12J 




3 


7 


49 


2' 


12' 


~20 


3 


7 


11 


2' 


12' 


30 



761 



(A.13) 



3. In NNLO we extract the feo-pi'oportional contribution by collecting all terms in the 
general structure of T2 computed in [3]. The result is 



T^'^(cu,x)+ln 



5! 



To ® {V. 



(l)^ 



ln(^^^T,(g^g^/i^.^)^^ ll^2^g^ [To®Vj(^,x)}, (A.14) 



where 



TP{u;,x)=To{Q-\q';y) 



T^'\y,x) + LNiu,y) (^j' V - T« iy, x) - ^LN\u,y)Viy,xy 



(A.15) 
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We use here for the elements of T2 the notation Tp\ Tft\ '^''^^ and Vn^\ introduced in [3], 



that differs from the previous one by a factor of 2. We have 



^(2) 



29 
12 



ya _^ y« 



209 



3 
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-9 + 9 



{x,y)-Q5{x-y), (A.16) 



while T^^^ is defined in (lA.Sbp . One observes that these T-kernels are calculated in terms 
of the kernels V and which enter the evolution kernel V^^\ and the derivatives V ^ V that 
will be discussed in Appendix |Bl The origin of g has been clarified in Sec. 3 of Ref . 0] • 

Under the assumption that T2 bo - Tp, it turns out that the entire LN(ix', y) dependence 
appears only inside the term Ln(?/) [cf. Eq. fl2.7p ]. in analogy to the NLO case, see flA.Tp : 



T2 



■ bo ■ T/s; 



(A.17a) 



In 



l4\ 



[r«(y,x) + Ln(y)y+(y,x)] 



- \ Ln^(2/)^(2/,a;)_ 



(A. 17b) 



This important property of the Tp structure has already been mentioned in 0]. While all 
terms in (lA.17bp contribute to the discontinuity of T^, only the last term, which contains 
the square of a logarithmic expression, contributes a new type of discontinuity. On the other 
hand, the term proportional to ln(yu|//i^) equals Ti, as it can be seen from Eq. (lA.7p . Thus, 
the final structure of Tp assumes the form 



+ CFTo(g^g^;y)®{T/)+Ln(y) [(\/«)+ - r«] - iLn2(y)K,}(y, x) . (A.18) 

From this expression, one sees that the first term, which contains no logarithm at all, has 
no influence on the discontinuity of Tq. On the other hand, the remaining terms contain 
logarithms, which do affect the discontinuity, with the last one being the new contribution 
first appearing at the NNLO level. 

Then, the partial amplitude Tp ® ipn reads 



Tp ® ^„ 



-In 



/4 



+CfTo ® i rj'^ + Ln 



-Ln^. (2^7(72)) (A.19) 



Let us now calculate the discontinuity entailed by the new contribution in NNLO, 

'{Q^ + s)y-s' 



IT 



1 



2{s + Q^) 7o y-yo 



s; y) (g) In^ 
1 dy 



{\n{y - yo) - \n{yo)) + In 



/^F 



(A.20) 
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Substituting identity jl] 



-Im 



TT 



ln^(?/ - 1/0 ) 



2/ -1/0 

into Eq. f[O0|) . we obtain for its RHS 
1 



Oivo > y) 





- 2/o) 


In^ 




-ttVS 












-{ 


^(z/o > y) 


\n{\y 








+(yo 


yo) 









in(i/ - yo) 
y-yo 



(A.21) 



(A.22) 



APPENDIX B: STRUCTURE OF THE EVOLUTION KERNEL V 

This Appendix compiles the crucial properties of the evolution kernels borrowing results 
from 65, 6^] and Following 67, 68|, we introduce the auxiliary kernels V'^''^{x, y; A), viz.. 



V:{x,y;X) 
Vl{x,y;X) 



9{y > x) 



l+A 



9{y > x) / X 
y-x \y 



l+A 



+ 



X — 


^ X 


y - 


^ y 


X - 


X 


y - 


y 



(B.l) 
(B.2) 



and their sum V{x,y; A) = V^{x,y; A) + V^{x,y; A) that includes the main of the logarith- 
mic contributions — generated by the one-loop renormalization of the running coupling — 
accompanied by the factor ln{x/y). Just the effects of this renormalization lead to those 
contributions that are proportional to Bq and pertain to the derivatives of the auxiliary kernel 
V{x,y;X), i.e.. 



2rix,y;X) \x=o= Vix,y) = 2 
2V"{x,y-X) U=o=V{x,y) = 2 



X 



eiy>x)-[l 

y \ y-x 



X 



eiy>x)-[l 

y \ y-x 



In 
In^ 





X - 


X 




[y- 


y 




X - 


X 




[y- 


y 



, (B.3) 
, (B.4) 



we already faced in Eq. flA.16|) . The ^o-proportional part Vg^^ of the NLO kernel V'^^\ that 



enters Tp in Eqs. (IA.17bp and flA.lSP for T^, contains the non-logarithmic terms 



Vr. 



(1) 



These contributions can be obtained from the generating kernel (66 



Vp{x, y\X) = [(1 + X)V^{x, y; X) + V\x, y; A)] C(A) 



(B.5) 



(B.6) 



for any order of Bq. Here C{X) is an analytic function in the variable A with C(0) = 1. 
To obtain the 6o-contribution at any desired fixed order of the parameter agbo, one has to 
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expand the kernel Vi3{x,y\asbo) in a Taylor series with respect to up to this order. Hence 
Vi3{x,y\0) = V+{x,y), while the first differentiation of Vg with respect to as, 



-^Vi3{x,y\asbo] 



as=0 



(B.7) 



leads to V^^^"*. The generalized kernel Va(a;, y\X) has been derived from the diagrams for the 
ordinary one-loop kernels, generated by replacing single gluon lines by a sum of renormalon- 
chain insertions. The coefficient C(A) in Eq. (1B.7P accumulates non-logarithmic parts of 
these renormalon-chain contributions and was determined in 



APPENDIX C: MAIN ELEMENTS OF THE NNLO PARTIAL AMPLITUDE 



We present now the main elements of the partial amplitudes ® ipn and 



ipn entering Eq. (lA.lQD in Appendix |Al We split each of the expressions 
Delow in two parts: the singular part with x ^ is extracted in an explicit form that 
turns out to be proportional to tpn, whereas the other part contains an integration over 
longitudinal momentum fractions. 

1. Recalling Eq. flA.Sbl) in App. \^ we find 



(r(i) ® ^n) (x) = 3 [1 + v\n)] Mx) + 



In^^U- 

X 



i^n{x) 



du 



y 



u — X 

3 + '"^'jj-T + 2 











f (x - 











X 



IT 



i^o{x) 



(C.lb) 



1 



Ln^ ■ V+ ® ip-n 



-Ln^ ■ v{n) ■ ipn 



--Ln^ ■ V"+ ® ^0 = 0. 



(C.2a) 
(C.2b) 



2. According to the definition of Vg in (IB. 51) and the definitions of V , in ( IB. 31) . we 
have 



V+^^n) (X) 



|r^(i+^)iu(i)i*.„ 

{Jx U \ U — X J \U/ 



(u) — 1pn{x)] du+ {x ^ x) 



X In^ (x) — xln^(x) -\ h In^ ( — ) 



7 



X\ TX 



6x ln(x) + 6x ln(x) + 



5 + In' 



X 



X\ TT 



X 



i^n{x) 

i'oix) 



(C.3a) 
(C.3b) 
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x] 



-v{n) +t;"(n) 



^Jnix) 



-3 



X ln(x) + X ln(x) 



5 + ln^ 



X 



X 



71 



(C.4a) 



Mx) , (C.4b) 



Jx U \UJ 



(u) - ilJn{x)]du + ilJn{x)x ( ^ - In^(x) 



+ [X ^ X) , 

® "^0 = 3 (x ln(x) + X ln(x) + 3xx) . 



Using the expressions obtained in Eqs. fIC.lbl) and flC.4bp . we obtain 

(2V^j^^)+ - r(^)j 0tpo = -6x\n{x) - 6xln{x) . 



(C.5a) 
(C.5b) 



(C.6) 



3. The most characteristic contribution to is represented by the term Tp'^\ The result 



for the convolution is 0] 



29 209 • 7 • 1 •■ 
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(C.7) 



This expression contains a couple of new kernel elements, notably, V+ and The convo- 
lution expressions for each of these kernels with the Gegenbauer harmonics are displayed 
below, starting with the general case n, and followed by the zeroth-order harmonic. 



VV ®^n) {X) 



u \ u — X 



X 



In ( - ) [ipn{u) - ipn{x)]du + {x ^ x) 



+ i^n{x) 



41n(x)Li2(x) -4Li3(x) xln^(x) + 21n(x) In^(x) 

3 



X 

+ 4 

2 



+ (x x) 



(C.8a) 



K|_ ® V^o ) (x) = 6xx|4 [ln(x)Li2(x) + ln(x)Li2(x) — Li3(x) — Li3(x)] 

2 

- -(In^(x) + ln^(x)) + 21n(x) ln(x)(ln(x) + ln(x)) + 6 
o 

— 6x (in^(x) + ln(x)) — Qx (in^(x) + ln(x)) , 



(C.8b) 



9+{x,y) 



djy - x) 
y-x 



X 

1 - - 

y 



L.(l)-i.n'(l-^ 



y^y. 



, (C.9) 
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2 [g+ ® i)n] (x) 











U — X 





Liofl 



X_, 

u' 

1 



Li2(l)--ln^ 1 

2 \ u 



du 



vr 



ln(x) In^(x) - - In^(x) - — ln(a;) + 4Li3(x) 
3 3 



- 6C(3) 



+ (x ^ x) 



(C.lOa) 



2 [qj^ ® t/'o] (2;) = 6xx| ln(x) In^(x) + In^(x) ln(x) — - [in'^(x) + In' 



X 



- — [ln(x) + ln(x)] + 4 [Li3(x) + Li3(x)] - 12C(3) - 2 
-12xln(x) - 12xln(x) . (C.lOb) 

Finally, substituting Eqs. (IXQ]) . (lOSbD . (laSbj) . flOSbl) . and (IC.lObD into Eq. flOTl) for n = 0, 

/ox 

and taking into account that i;"(0) = f (0) = 0, one obtains for T^^ ® ipQ the expression 
rj2) = xxjsO [Li3(x) + Li3(x)] - 6 [ln(x)Li2(x) + ln(x)Li2(x)] - [in^(x) + ln^(a 
-5 In^ (^-) + [ln(x) + ln(x)] (3 ln(x) ln(x) - 27r^) - 72C(3) - -tt^ - ?} 



19 3 
+ — [xln(x) + xln(x)] + - [xln^(x) + xln^(x 
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